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AIRS/AMSU on Aqua 

AIRS/AMSU launched on EOS Aqua May 5, 2002 

AIRS is a multi-detector array grating spectrometer 

2378 channels between 650 cm-1 and 2760 cm 1 

Channel spacing = v/2400 (0.25 cm 1 - 1.1 cm' 1 ) 

Resolving power v/Av = 1200 (0.5 cm 1 - 2.2 cm' 1 ) 

AIRS Field of View (FOV) 13 km x 13 km at nadir 

AMSU is a microwave sounder with a nadir 45 km x 45 km Field of 
Regard (FOR) 

AMSU-A has flown on previous satellites 
3x3 AIRS FOV's within an AMSU-A FOR 
One AIRS/AMSU sounding is produced per AMSU-A FOR 
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Overview of AIRS/AMSU Retrieval Methodology 

Physically based retrieval system 

Independent of GCM except for surface pressure - used to compute 
expected radiances 

/\ 

Uses cloud cleared radiances R| to sequentially determine the solution for 

different geophysical parameters 

/\ 

Ri represents what AIRS would have seen in the absence of clouds 
All AIRS/AMSU retrievals have error estimates which are used for Quality 
Control 

Goddard DISC has been analyzing AIRS/AMSU data using AIRS Version-5 
algorithm 

Retrievals are near real time 

Analyzed data from September 2002 through the present 
AIRS Version-6 algorithm is finalized - see paper 8510-29 Wed. at 11:20 a.m. 
AIRS Version-6 will become operational in late 2012 for forward processing 
All previous data will be re-analyzed to produce a consistent climate 
data set 
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AIRS/AMSU Products 

Surface and Atmospheric Products - one set per AMSU FOR (45 km) 

Sea/land skin temperature T s 
Temperature profile T(p) to 1 mb 
Water vapor profile q(p) to 100 mb 
0 3 , CO, CH 4 ,C0 2 profiles 
Cloud Cleared Radiances Rj 

Cloud Products - one set per AIRS FOV (13 km). 

Radiatively effective cloud fraction ae for up to two cloud layers 
a is geometric fractional cloud cover 
s is cloud emissivity at 11 pm 
Cloud top pressure p c for up to two cloud layers. 

Outgoing Longwave Radiation (OLR) is computed using ae, p c , and 
retrieved parameters. 

Cloud clearing theory says that 

1) Longwave channels should be used to determine Rj for all channels. 

2) T(p) and T s should be determined using R; only in shortwave 

channels and longwave channels not sensitive to clouds. 
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Sample AIRS Cloud Free Brightness Temperature 
AIRS Channels Used in Different Retrieval Steps 


(K) 

Woven umb e r 
650 to 785 




cm 


( v \ 2 90 

( ^ 275 

Wavenumber 260 

785 to 1150 245 

230 

215 

3ol* 

290 
275 

Woven umb e r 2 ®^ 

1200 to 1650 245 

230 

215 

30h 2m 

290 
275 
260 


cm 


(K) 


cm 


(K) 

Woven umb e r 
2175 to 2700 
cm -1 


245 



* C I o u d 

* Wa ter 

* CH 4 


Clearing 

Vapor 


* T emp erature Profile *C0 

* Ozone 

* LW Em i s s i v i t y 


National Aeronautics and Space Administration 


Joel Susskind, Gyula Molnar, Lena Iredell, Robert Rosenberg , Oreste Reale - SPIE Paper #851 1-25 


5 


Objectives of AIRS/AMSU 


• Provide real time observations to improve numerical weather prediction. 

Could be R, (used by NCEP, ECMWF) orT(p), q(p). 

/\ 

Accuracy of R; T(p), q(p) degrades slowly with increasing cloud fraction. 
There is a trade-off between accuracy and spatial coverage. 

Using soundings or radiances only in clear cases limits utility of the 
data. 

• Provide observations to measure and explain interannual variability and 
trends. 

Must provide very good spatial coverage but also be unbiased 
Can be less accurate than needed for data assimilation. 

• Error estimates and quality flags provide options for use in either 
weather or climate applications. 
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AIRS/AMSU Version-5 Quality Control 


All AIRS/AMSU retrievals and cloud cleared radiances have level-by-level and 
channel-by-channels error estimates. Thresholds of error estimates are 
used for Quality Control (QC) of all retrievals - determines pressure p best 
down to which T(p) is acceptable. 

Version-5 used Standard QC Thresholds, optimized to give a best middle 
ground for Data Assimilation (highest accuracy) and Climate (best spatial 
coverage). Tighter QC thresholds perform better for Data Assimilation (DA). 

Much outstanding science has been conducted using AIRS Version-5 
products. 

I will show two studies related to: 1) improving forecast skill and 
2) explaining why global and tropical mean OLR has decreased recently. 

Version-6 will have further improvements, including separate sets of 
thresholds optimal for weather and climate applications and an 
"AIRS Only" (AO) mode to serve as a back-up should AMSU-A continue to 
degrade. 
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PRESSURE (mb) 


Global Temperature Profile Statistics 
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Version-6 T(p) retrievals with DA QC have RMS errors < IK throughout troposphere 

Version-6 T(p) retrievals with Climate QC have much greater yield than Verison-5 with small biases 

Differences between Version-6 and Version-6 AO are small 
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Data Assimilation Experiments 

Ran experiments with GMAO GEOS-5 Data Assimilation System (DAS). 

• Uses NCEP Operational GSI Data Assimilation Methodology. 

• Forecasts and NCEP GSI Analysis run at 0.5° x 0.625° resolution. 

We compared forecasts from four assimilations: 

Control No AIRS 

AIRS V5 temperature profile, Standard QC 

AIRS V5 temperature profile, Tight QC 

AIRS radiance assimilation (Operational) 

Control used all the data NCEP used operationally in 2003. 

• Assimilates all satellite data but AIRS, including Aqua AMSU radiances 

AIRS Radiance Assimilation (Operational procedure) adds observed airs 
radiances unaffected by clouds to Control - excludes most AIRS observations. 

Control + AIRS (Standard or Tight) adds Version-5 global Quality Controlled 
temperature profile retrievals to Control. 

• Assimilated as if radiosonde data. 

• Error estimate of temperature profile is used as the level measurement error. 
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Forecast Impact Tests 


Ran four sets of experiments, covering different seasons and years. 

January 1 - January 31, 2003 

October 15 - November 19, 2005 

August 10 - September 16, 2006 

April 15 - May 18, 2008 (Tropical Storm Nargis) 


The AIRS temperature profiles are assimilated down to p best determined 
using appropriate QC thresholds. 

Seven day Forecasts run from each 0 Z Analysis for each experiment. 

The accuracy is judged against anomaly correlation of 7-day forecasts 
vs. ECMWF Analysis for that time. 
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Anomaly Correlation 


500mb Geopotential Heights 

Northern Hemisphere Extra Tropics 



GEOS-5 Control (No AIRS) 

G EOS-5 AIRS Standard 

GEOS-5 AIRS Tight 

GEOS-5 Radiance (Operational) 


Day 

Average of Summer 2006, Spring 2008, Fall 
2005, and Winter 2003 


An anomaly correlation of 1.0 represents a perfect forecast 

An anomaly correlation of 0.6 is the lower bound of a useful forecast 

AIRS Tight improves seven day forecast skill by about five hours compared to No AIRS 
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Example of AIRS Impact on Forecast Track 


Tropical storm Nargis 
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Analysis 
Well-defined 
Cyclone 

Green: 

Observed 

Track 

7SE 80E 82E t IE 86E 88E 90E 92E 94E 95E 98E "7SE 80E 62E 84E 86E 98E SOE ME 9<E 96E SSBE 


CNTRL 06z 28April (6h FC) AIRS 06z 28April (6h 




CNTRL 12z 2 May ( 1 0Sh FC) 



AIRS 1 2z 2May (1 08h FC) 



AIRS 108-hour 
Forecast (sip) 


Green: 

Observed 

Track 


CNTRL Analysis (above) Accurate landfall is produced in the forecasts initialized with 

And forecast (below): No Cyclone AIRS: (Reale et al., 2009, Geophys. Res. Lett.) 
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7-day accumulated precipitation forecast initialized at July 22, 2010 OOz 

GEOS-5 Forecast using v5 Clear-sky radiances (RAD) GEOS-5 Forecast using v5 Retrievals (RET) 
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Assimilation of AIRSv5 retrievals goes in the direction of the observations 
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First Experiment Using NCEP GFS 

Lou Uccellini, Director of NCEP, was very impressed with the potential of 
this methodology to improve operational forecast skill. 

New methodology must be first tested with the GFS system used 
operationally at NCEP. 


Lou Uccellini, together with Tsengdar Lee (HQ) arranged to have the NCEP 
GFS ported to GSFC. 

We have just started to run the first experiment assimilating AIRS 
Version-5 T(p) in GFS. 

Covers the period July-August 2010, containing Pakistan flood. 

Results are encouraging but too preliminary to be shown at this time. 

Future experiments will contain 

• Tests with other time periods. 

• Optimization of DA QC thresholds for use in GFS. 

• Experiments assimilating Version-6 T(p) retrievals. 
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Example of Climate Research 

• Comparison of AIRS and CERES anomaly time series of OLR. 

• Explanation of recent decreases in global and tropical mean values of 
OLR. 

Data sets used 

AIRS Science Team Version-5 monthly mean data obtained from 
Goddard DISC (level-3). 

OLR, q 500 , cloud fraction ot£. 

Data products extend to July 2012. 

CERES Science Team monthly mean data obtained from Langley ASDC 

CERES Terra EBAF (Energy Balanced And Filled) Edition 2.6r. 

Data products extend to June 2011. 

CERES is considered the "gold standard" for OLR 

AIRS/CERES comparisons are shown for the period September 2002 
through June 2011. 
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Definition of Anomalies, ARCs and ENCs 


Eight-year monthly climatologies were generated for each l°x 1° grid box 
by averaging data for eight Januaries, eight Februaries, 

The monthly anomaly for each grid box is the difference of the monthly 
mean value for that month from its climatology. 

The Average Rate of Change (ARC) for a grid box is the slope of the 
straight line passing through the monthly anomaly time series. 

Values of ARCs depend on the extent of the time series used. 

Spatial patterns are more important than precise values. 

An area mean ARC is the cosine latitude weighted average ARC over the 
area. 

The El Nino Correlation (ENC) for a given grid point is the temporal 
correlation of the anomaly time series for that grid point with the 
El Nino Index. ENCs should be less time period dependent. 
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OLR Anomaly Time Series 
September 2002 through June 2011 


Global Mean 


Tropical Mean 




CERES 

AIRS minus CERES 


CERES 

AIRS minus CERES 

1.5 times the El Nino Index 

El Nino Index ARC = -0.123 ± 0.046 K/yr 
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OLR Anomaly Time Series Comparison 
September 2002 through June 2011 


Data Set 

Global 

Tropical 

AIRS ARC (W/m 2 /yr) 
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AIRS and CERES global mean and tropical mean anomaly time series agree very closely 
with each other. 

Both show that global, and especially tropical, mean OLR have decreased over the time 
period under study. 
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Important Highlighted Regions in Subsequent Figures 


NOAA Nino-4 Region: 

We define the El Nino Index as the 
NOAA monthly mean T skin anomaly, 
compared to 8-year climatology, averaged 
over the NOAA Nino-4 region 
5°N to 5°S, 150°W, westward to 160°E 




, VSr- 


Warm Pool Maritime Continent (WPMC) region: 

West of Nino-4 region 


Equatorial Eastern Pacific and Atlantic (EEPA) region: 

East of Nino-4 region 


National Aeronautics and Space Administration 


Joel Susskind, Gyula Molnar, Lena Iredell, Robert Rosenberg , Oreste Reale - SPIE Paper #851 1-25 


19 




ARCs and ENCs of OLR Anomaly Time Series 
September 2002 through June 2011 
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b) CERES ARCS 
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b) CERES ENCs 



m i i i ni 

-1.0 -.84 -.68 -.52 -.36 -.2 -.04 .04 .2 .36 .52 .68 .84 1. 

Global Mean=0.03 STD=0.36 


OLR ARCs are very negative over the WPMC and EEPA regions (blue, yellow). 
OLR ENCs are very positive over the WPMC and EEPA regions (blue, yellow). 
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Area Mean Statistics for AIRS and CERES OLR 
September 2002 through June 2011 

AIRS CERES 


Spatial Area 

OLR ARC 
(W/m 2 /yr) 

OLR ENC 

OLR ARC 
(W/m 2 /yr) 

OLR ENC 

WPMC Region 

-1.502 ±0.525 

0.871 

-1.496 ±0.529 

0.870 

Tropical outside 
WPMC Region 

0.004 ± 0.054 

0.050 

0.034 ± 0.053 

0.101 

Global outside 
WPMC Region 

-0.030 ± 0.028 

-0.129 

0.005 ± 0.027 

-0.372 

EEPA Region 

-0.631 ±0.158 

0.767 

-0.611 ±0.154 

0.761 

Tropical outside 
EEPA Region 

-0.037 ± 0.048 

0.599 

-0.011 ±0.047 

0.511 

Global outside 
EEPA Region 

-0.044 ± 0.020 

0.256 

-0.011 ±0.019 

0.039 


ARCs of otherwise Tropical mean and Global mean OLR computed 
outside of either WPMC region or EEPA region are essentially zero. 


National Aeronautics and Space Administration 


Joel Susskind, Gyula Molnar, Lena Iredell, Robert Rosenberg , Oreste Reale - SPIE Paper #851 1-25 


21 




































AIRS Products September 2002 through June 2011 
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T skin ARCs and ENCs are for the most part of opposite sign (similar color). This shows most ARCs 
are El Nino driven, except near the North Pole. 

ENCs of q 500 and as are similar to each other and are both very negative over the WPMC and 
EEPA regions. This results in positive OLR ENC's in these regions. 
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Findings of Climate Study 

Recent AIRS and CERES OLR anomaly time series are in very close 
agreement on a 1° spatial scale. 

Both show a significant decrease in global mean and tropical mean OLR 
over the period September 2002 through June 2011. 

The degree of agreement is remarkable, but the recent OLR decrease is not 
because it is the result of a transition from El Nino conditions to La Nina 
conditions. This implies nothing about long term changes. 

The recent decrease in global and tropical mean OLR results from the very 
large negative correlations of water vapor and cloud cover in each of the 
WPMC and EEPA regions with the El Nino Index. This more than 
compensates for the reverse effect in the Nino-4 region. 

AIRS is expected to last through the life of Aqua =2022. 

This will allow for a 20 year AIRS Climate Data Set. 

Version-6 climate products will be further improved over Version-5. 
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